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ABSTRACT 

Equations  are  developed  to  estimate  velocity  turbulence  from  index  of  refraction  turbulence  as  measured  by 
the  thermosonde.  Using  an  adaptation  of  the  radar  “power”  method  to  estimate  velocity  fluctuations  from 
backscattered  radar  return  and  concurrently  measured  atmospheric  data,  the  method  is  applied  to  optical 
turbulence  measured  by  the  thermosonde.  The  method  is  then  used  to  compare  velocity  turbulence  levels 
as  estimated  by  the  thermosonde  to  those  estimated  by  the  “GroundWinds”  lidar  in  the  lee  of  Mt. 
Washington  in  New  Hampshire.  There  is  reasonable  agreement  in  identification  of  the  areas  of  higher 
turbulence. 


NOMENCLATURE 
A  =  Constant  or  coefficient 
=  Constant  or  coefficient 
^  =  Structure  constant  of  the  quantity  x 

[units  of 

(r)  =  Structure  function  of  quantity  x  [units  of 

X  squared] 

=  Distance  [m] 

=  Acceleration  of  gravity  [m/s^] 

=  Ratio  of  specific  heats  Cp/Cy  [n.d.] 

=  Spatial  wave  number  [m"^] 

=  Gradient  of  refractive  index  [m“^] 

=  Brunt-Vaisala  Frequency  [s'^j 
=  Index  of  refraction  [n.  d.] 

=  Pressure  or  partial  pressure  when 
subscripted  [Pa] 

=  Specific  humidity  [n.d.] 

[k)  =  One-sided,  one-dimensional  spatial  power 

spectrum  of  x 
=  Temperature  [K] 

V  =  Components  of  mean  horizontal  wind 
velocity  [m/s] 

=  Fluctuating  part  of  horizontal  wind  [m/s] 

=  Distance  (m) 

=  Altitude  in  Plots  (km) 
z  =  Distance  in  vertical  direction  (m) 

e  =  Turbulent  kinetic  energy  dissipation  rate 

[Wkg-'ormV^] 


y  =  Constant  in  Cn^  -  e  relationship  [n.d.] 

X  =  Wavelength  of  refracted  radiation  [m] 

0  =  Potential  temperature  [K] 


Subscripts: 

d 


n 

T 


U  or  u 
wv  = 
Superscripts: 
Bar 


Dry  air 

Index  of  refraction 
Temperature 
Velocity 
Water  vapor 

Time  averaged  or  mean  component 


INTRODUCTION 

Turbulence  in  the  atmosphere  is  difficult  to  forecast  and 
difficult  to  measure.  Atmospheric  turbulence  is  the 
random  fluctuation  of  the  air  mass  caused  by  shear 
instabilities  in  jet  streams  or  over  mountains,  by 
convective  instabilities  associated  with  thunderstorms, 
and  by  the  breaking  of  buoyancy  or  “gravity”  waves, 
especially  at  higher  altitudes.  The  spatial  structure  of 
atmospheric  turbulence  is  often  characterized  by  the 
rate  of  turbulent  kinetic  energy  dissipation  per  unit 
mass,  often  called  the  eddy  dissipation  rate,  e.  While  it 
refers  to  the  conversion  of  turbulent  kinetic  energy  to 
heat  at  the  molecular  level,  it  is  implicitly  tied  to  the 
magnitude  of  fluctuations  throughout  the  range  of 
turbulent  scale  sizes  through  turbulence  theory. 
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Clear  Air  Turbulence  (CAT),  the  effect  of  strong  waves 
and  turbulence  in  the  atmosphere  on  air  vehicles,  is  one 
particularly  dangerous  manifestation  of  turbulence. 
Forecast  of  CAT  is  the  subject  of  continuing  research 
by  NOAA,  NCAR  and  the  U.S.  Air  Force  and  Navy. 
Many  detection  algorithms  have  been  evaluated,  and 
generally  each  algorithm  exhibits  a  wide  variation  in 
performance.  This  is  not  only  due  to  the  complexity  of 
the  forecast  problem,  but  also  the  means  available  to 
validate  the  performance,  primarily  pilot  reports 
(PIREPs).  These  PIREPs  depend  on  an  aircraft  being  at 
the  right  place  at  the  right  time.  The  perceived  severity 
of  a  CAT  event  is  very  dependent  on  the  specific 
aerodynamic  characteristics  of  the  aircraft  and  the 
subjective  judgment  of  the  pilot.  Also  the  events  that 
usually  get  reported  are  at  Ihe  more  severe  end  of  the 
possible  spectrum  of  turbulence. 

TURBULENT  KINETIC  ENERGY 

At  equilibrium,  Kolmogorov  hypothesized  that 
turbulent  kinetic  energy  cascades  from  the  largest 
eddies,  called  the  outer  scale  of  the  turbulence,  to  the 
smallest  eddies,  the  inner  scale,  primarily  by  inertial 
processes\  At  the  inner  scale,  kinetic  energy  is 
converted  to  thermal  energy  by  viscous  dissipation  at 
the  dissipation  rate,  S .  By  virtue  of  the  Kolmogorov 
hypothesis,  small-scale  isotropy  is  also  implied.  When 
the  velocity  spectrum  is  within  the  inertial  range 
(between  the  outer  and  inner  scale),  Kolmogorov 
reasoned  by  dimensional  arguments  that  the  one  way, 
one  dimensional  power  spectral  density  of  longitudinal 
velocity  fluctuations  has  the  following  functional 
relationship^: 

(1) 

where  A  is  a  constant  equal  to  approximately  V2.  In  the 
inertial  range,  the  second  order  cross-stream  and 
vertical  components  of  velocity  structure  functions  in 
the  longitudinal  direction  are  larger^  by  a  factor  of  4/3, 
In  the  large  Reynolds  numbers  associated  with  the 
atmosphere,  there  can  be  a  wide  inertial  range  over 
which  Equation  1  applies.  Turbulent  equilibrium 
implies  that  a  constant  level  of  turbulent  kinetic  energy 
has  been  reached  and  stays  steady,  £“  is  a  constant,  and 
the  power  spectral  density  stays  constant  within  the 
inertial  range  and  appears  as  a  straight  line  of  slope  -5/3 
on  a  log-log  plot.  The  variance  of  the  turbulent  energy 
fluctuations  in  longitudinal  velocity  within  the  inertial 
range  is  determined  by  the  equation: 


(2) 

The  spectra  can  also  be  expressed  in  terms  of  the 
velocity  structure  constant, 

S,Xk)^0.25Clk-^'^  (3) 

By  comparison  of  Equations  1  and  3,  , 

where  a  is  four  time  A,  or  approximately  2. 

The  radar  provides  a  means  of  remotely  estimating 
turbulence  on  a  continuous  basis.  There  are  two  basic 
methods  of  estimating  the  velocity  turbulence  using 
radars'^,  one  from  reflected  power  and  the  other  using 
the  spectral  spread  of  the  returned  beam.  The 
methodology  of  the  first  technique  has  been  adapted  to 
the  thermosonde.  In  this  paper  we  briefly  introduce  the 
radar  backscatter  or  power  method  of  obtaining 
turbulence  strength,  then  describe  the  thermosonde 
instrument  and  the  adaptation  of  the  power  method  to 
the  thermosonde.  Finally  estimates  of  the  velocity 
turbulence  are  presented  from  a  campaign  in  Bartlett, 
New  Hampshire  and  compared  to  turbulence  levels 
estimated  using  the  Ground  Winds  Lidar  in  Bartlett. 

THE  RADAR  POWER  METHOD 

The  index  of  refraction,  n,  in  the  atmosphere  can  be 

written  as^: 

n  =  \+aj{Z)PJT+a^(A)P^/T  (4) 

where  P  is  the  partial  pressure  of  either  the  dry  air  (d) 
or  the  water  vapor  (vw)  and  the  a  coefficients  are  the 
wavelength  dependent  coefficients  for  either  the  dry  air 
or  the  water  vapor.  Fluctuations  in  the  index  of 
refraction  are  caused  by  fluctuations  of  either  the 
density  of  the  dry  air,  the  concentration  of  water  vapor 
in  the  air,  or  both.  Radars  detect  fluctuations  in  the 
index  of  refraction  of  the  view  volume  by  Bragg 
scatter^.  For  VHF  radars  the  backscatter  will  generally 
fall  within  the  inertial  sub  range  of  the  turbulence. 

The  basis  of  the  reflected  power  method  is  that  the 
fraction  of  transmitted  power  that  is  reflected  back  to 
the  radar  is  proportional  to  the  mean  refractive  index 
structure  constant,  Therefore,  the  method  can  only 
be  used  for  calibrated  radars.  The  measured  C„^  is  the 
mean  value  over  the  sampled  volume  of  the 
atmosphere,  which  is  determined  by  the  beam  spread 
and  the  range  gates  of  the  radar.  Estimates  of  the 
maximum  Cn^  within  the  volume  is  then  based  on  an 
estimate  of  the  fraction  of  the  volume  that  is  actually 
turbulent.  The  actual  form  of  the  radar  equation 
depends  on  the  radar  configuration"^. 
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The  relationship  of  Cn  to  s  invokes  equilibrium  and 
the  relationships  for  the  dissipation  of  turbulent  kinetic 
energy  and  refractive  index  turbulence^’  Hocking 
and  Mu^  consolidated  the  relevant  considerations  into 
the  formula: 


,2  _  1 


(5) 


where  N  is  the  Brunt  Vaisala  frequency  and  M  (at  radar 
wavelengths)  is  the  refraction  index  gradient,  which  are 
defined  by  the  following  equations: 


'  e  dz 


(6) 


M^-TT.bxlO-" 


1  I5500q 

1  dh.qldz^ 

\  &  J 

^  T 

^  2d\n.6ldz) 

(7) 


where  6  is  the  potential  temperature  defined  by  the 
equation^: 


^(Z)  =  r(Z) 


P(Z)^ 


(8) 


Ps  is  the  reference  pressure,  usually  1000  hPa.  Finally, 
Y  in  Equation  5  was  assumed  constant  in  the  early 
literature,  but  is  often  expressed  as  a  function  of 
gradient  Richardson  Number^,  defined  by: 


where  U  and  V  are  mean  horizontal  wind  components. 
Some  of  uncertainties  of  the  radar  power  method  arise 
from  the  uncertainties  in  the  calibration  of  the 
microwave  source,  the  estimate  of  the  fraction  of  the 
atmosphere  that  is  turbulent,  the  form  of  the  / 
parameter,  and  the  requirement  to  rely  on  available 
balloon  soundings  to  determine  the  temperature  and 
humidity  profiles,  the  gradient  Richardson  Number  and 
the  Brunt“Vaisala  frequency^®.  Using  balloon-borne  in 
situ  temperature  and  velocity  turbulence  measurements, 
Bertin,  et  al*^  deduce  that  the  beam  spread  method, 
while  suffering  from  some  uncertainties,  is  probably 
more  accurate  than  the  power  method  that  is  used  here. 
There  is  an  extensive  database  of  temperature  variance 
measurements  using  the  thermosonde.  These  data  have 
also  been  used  to  estimate  the  velocity  turbulence^’ 


THERMOSONDE  MEASUREMENTS: 

The  thermosonde  is  designed  to  measure  optical 
turbulence  in  the  atmosphere,  that  is  refraction 
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turbulence  in  the  visible  or  near  visible  wavelengths. 
For  these  wavelengths,  the  magnitude  of  the  water 
vapor  term  in  Equation  4  is  so  small  that  it  hardly 
contributes  to  the  index  of  refraction,  except,  perhaps, 
at  very  low  altitudes  over  a  body  of  water.  While  both 
pressure  and  temperature  appear  in  Equation  4, 
fluctuations  in  pressure  in  the  free  atmosphere  are  small 
and  they  dissipate  rapidly  through  acoustic  processes. 
Therefore,  in  situ  measurement  of  optical  turbulence  is, 
in  practice,  reduced  to  measurement  of  fluctuations  in 
temperature. 

Before  covering  the  calculation  of  e,  we  give  a  brief 
description  of  the  thermosonde  as  implemented  at  the 
Air  Force  Research  Laboratory*^’^"^.  Two  wire  probes, 
spaced  one  meter  apart  on  a  horizontal  boom,  sense 
temperature  differences  as  the  instrument  ascends 
through  the  atmosphere  at  around  5m/s,  towed  by  a 
balloon  on  a  1 10m  line.  The  probes  use  3.45iim 
diameter  wire,  which  has  a  time  constant  of  less  than 
1ms,  at  nominal  ascent  speeds  and  sea  level  conditions. 
The  temperature  difference  is  continuously  averaged  by 
an  onboard  analog  RMS  integrated  circuit.  The  output 
of  the  RMS  chip  is  transmitted  back  to  the  groimd 
station  at  1 .2s  intervals.  Since  the  balloon-borne 
payload  is  attached  to  a  meteorological  radiosonde, 
concurrent  wind  velocity,  pressure,  mean  temperature, 
and  humidity  information  are  also  relayed  to  the  ground 
station.  The  output  of  the  RMS  chip,  when  converted 
to  the  RMS  temperature  difference,  is  a  running 
structure  function.  Assuming  that  the  turbulence 
follows  the  Kolmogorov  hypothesis,  the  structure 
function  has  the  following  form^^: 

D^{d)  =  Cld'^'^  (10) 

where  d  is  the  distance  between  the  points  of 
measurement  and  Dj{d)  is  the  structure  function  and 
Ct  is  the  temperature  structure  constant.  Since  the 
thermosonde  acquires  temperature  differences  at  a  Im 
separation,  the  magnitude  of  the  Im-structure  function, 
Z)7’(lm),  is  the  magnitude  of  the  structure  constant.  The 
amplitude  of  balloon  acquired  temperature  structure 
constant  (C/)  typically  varies  from  the  thermosonde 
noise  floor  (1x10"^  to  IxlO'^K^m^'"^)  to  a  maximum  of 
around  O.OIKWI  Values  below  IxlO'^K^/m^^^do  not 
significantly  affect  optical  propagation.  Conversion  of 
the  temperature  structure  constant  to  the  refractive 
index  structure  constant,  C„^,  depends  on  local  pressure 
and  temperature  and  the  wavelength  of  the  radiation 
that  is  being  propagated.  For  radiation  near  the  visible 
spectrum,  when  the  contribution  can  safely  be  ignored 
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except,  perhaps  right  over  a  body  of  water,  it  is 
customary  to  use  the  formula*^*: 

Cl=(aM)PIT^)'Cl  (11) 

where  ad(X)=79xlO‘^KyPa  for  visible  and  near  IR 
wavelengths. 

The  computation  of  ^from  the  thermosonde  eliminates 
some  of  the  shortcomings  of  the  comparable  radar 
calculation.  First,  the  response  of  the  instrument  is  fast 
enough  to  eliminate  the  need  to  estimate  the  turbulent 
fraction  of  the  atmosphere  (it  is  unity).  Second,  the 
atmospheric  sounding  data  are  taken  concurrently  with 
the  Ct  data.  In  addition,  the  calculation  for  the 
thennosonde  does  not  require  terms  that  include  the 
contribution  of  water  vapor  to  the  index  of  refraction 
for  radar  wavelengths.  Rather  than  computing  8 
directly  from  the  temperature  structure  constant  (as 
shown  in  Ref  12)  we  will  simply  extend  the  radar 
power  method  (as  described  in  Ref. 8).  From  Equation 
5: 


values  of  Ri.  We  have  sided  with  other  authors^’^’^^ 
who  have  argued  that  gamma  approaches  a  non-zero 
constant  at  the  higher  values  of  Ri,  and  we  have 
maintained  a  value  of  y  equal  to  1  above  —  0. 12. 


Figure  1:  Plot  of  the  Gamma  factor  versus  the 
gradient  Richardson  number.  As  Ri  decreases,  the 
increasing  values  from  Equation  15  are  used.  Above 
0.12,  y  is  set  to  one. 


U  ’<  J  (12) 

For  near  visible  wavelengths,  humidity  can  be  ignored, 
and  M  simplifies  to: 


M  =  -aAX)^^ 


P_^ 

re  dz 


(13) 


At  this  point,  a  discussion  of  the  constant,  y,  is  in  order. 
Ottersten^  suggested  the  following  form  for  y: 


1  /"Pr-i?/"! 

^“a'Prl  Ri  j 


(14) 


Experimental  data  suggest  that  for  Ri«l,  the  turbulent 
Prandtl  number,  Pr,  is  near  unity.  We  follow  Hocking 
and  Mu®  in  keeping  Pr =1  and  use  their  expression  for 

y: 


^  a\  Ri  J 


(15) 


We  also  use  H&M’s  value  for  the  constant  ~  7.33. 
When  they  used  this  equation  to  analyze  thermosonde 
data,  they  used  Equation  1 5  for  values  of  Ri  less  than  1 . 
They  maintained  y  at  zero  for  higher  values  of  Ri.  This, 
they  argued,  was  to  suppress  regions  of  decaying 
turbulence,  which  should  be  the  case  for  the  higher 


*  In  Reference  12,  one  of  the  square  symbols  (^)  was 
inadvertently  omitted.  Please  correct  it! 
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As  a  final  observation,  when  the  equations  for  Cn  ,  N  , 
and  M  are  substituted  into  Equation  12,  we  obtain  the 
working  equation: 

-|3/2 

(16) 

dz 

Ultimately,  the  vertical  variability  in  the  estimated  s  is 
due  to  the  variability  in  both  Cj  and  the  derivative  of  0, 
with  the  former  varying  at  the  highest  frequency. 

THE  EXPERIMENT 

On  31  May  2001  thermosonde  instruments  were  flown 
adjacent  to  the  GroundWinds  Lidar  facility  at  the 
Bartlett,  NH.  This  site  is  24km  south  and  slightly  east 
of  the  summit  of  Mt.  Washington,  which  is  1,917m 
above  mean  sea  level  (MSL).  The  GroundWinds  site  is 
166m  MSL.  The  experiment  began  at  dusk,  ~9pm 
EDT,  with  two  thermosondes  launched  within  1  minute 
of  each  other.  The  next  two  thermosondes  were 
launched  right  after  1 1pm.  The  best  combination  of 
lidar  and  thermosonde  data  was  the  4*^  launch  of  the 
evening.  A  composite  plot  of  mean  air  temperature  and 
relative  humidity,  optical  turbulence  (C/),  wind  speed 
and  direction  for  the  flight  is  shown  in  Figure  2.  The 
results  are  quite  distinct  from  “typicaP’  thermosonde 
results,  but  they  may  be  representative  for  flow  in  the 
lee  of  a  good  size  mountain.  In  particular,  up  to  about 
4km,  the  optical  turbulence  levels  are  well  above  the 
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Clearl  model,  the  smooth  curve  in  the  Cn  plot.  At 
4km,  twice  the  height  of  the  mountain,  the  level  drops 
dramatically  to  well  below  Clearl  until  about  10km, 
perhaps  2km  below  the  altitude  of  minimum 
temperature,  which  is  the  start  of  the  tropopause  region. 

The  lower  4km  is  quite  active,  and  a  good  location  to 
expect  higher  turbulence  levels.  An  expanded  view  of 
the  lowest  5km  of  the  composite  plot  is  shown  in  Figure 
3,  The  horizontal  lines  in  the  wind  direction  plot  are 
“wrap  around”  lines  as  the  wind  direction  changes 
between  slightly  less  and  slightly  greater  than  true 
North,  and  not  an  indication  of  large  wind  sweeps.  The 
largest  region  of  turbulence,  above  the  expected  high 
activity  in  the  immediate  vicinity  of  the  ground,  is 
between  1 .5  and  2.0km,  near  the  altitude  of  the  top  of 
Mt.  Washington  (1 .9km).  This  could  be  the  result  of  a 
Kelvin-Helmholtz  instability  due  to  the  shear  over  the 
mountaintop.  The  high  turbulence  region  is  also  right 
above  a  local  maximum  in  humidity  at  about  70%, 
which  is  the  start  of  an  inversion  in  the  temperature 
profile  that  persists  through  this  Vi  km  band.  It  is  not 


Figure  4.  Weather  analysis^^  showing  a  jet  stream 
passing  through  northern  New  England  at  the  time 
of  the  test.  The  300hPa  altitude  contours  show 
approximate  jet  streamlines,  and  the  velocity  vectors 
show  relative  speed  and  direction  at  summit  height. 

unusual  to  have  increased  turbulence  at  the  top  of  a 
moisture  layer.  From  2km  to  3km  the  level  is 
somewhat  reduced,  followed  by  a  region  of  increased 
Cn^  between  3  and  V/t  km.  This  higher  region  is 
located  the  start  of  a  very  dry  layer.  Then  at  4km,  the 

5 


Cn^  level  drops  by  3  orders  of  magnitude  and  enters  the 
quiet  region  of  the  troposphere  discussed  above. 

The  quieter  region  is  actually  in  the  lower  part  of  the  jet 
stream,  shown  in  Figure  4  depicting  the  winds  over 
New  England  at  that  time.  Referring  to  the  wind  speed 
part  of  Figure  2,  the  jet  starts  with  a  high  shear  region 
between  3km  and  above  4km,  nearly  coincident  with 
the  second  region  of  high  optical  turbulence.  The  jet 
then  increases  to  nearly  40m/s  at  around  10km,  then  it 
decreases  to  about  20m/s  at  16km,  and  then  decreases 
rapidly  to  about  5m/s.  There  is  another  region  of  higher 
turbulence  around  the  upper  jet  shear  region  near  15km. 

RESULTS 

The  Cn  data  were  then  converted  into  estimates  for  6 
following  the  equations  presented  above.  The  full 
profile  of  dissipation,  and  the  intermediate  steps  of 
and  Ri  are  shown  in  Figure  5.  The  large  oscillations  in 
are  solely  due  to  the  small-scale  changes  in  the 
vertical  slope  of  the  potential  temperature,  0,  which  in 
turn  is  due  to  the  changes  in  slope  of  the  temperature, 

T.  Ri  is  then  obtained  by  dividing  by  the  square  of 
the  vertical  wind  shear.  The  generally  high  Gradient 
Richardson  Number  indicates  a  generally  stable 
atmosphere,  and  results  in  the  y  function  being  held  at 
unity  for  most  of  the  flight.  Only  in  the  lower  altitude 
regions  below  3km,  and  spikes  at  around  7  and  10km 
does  Ri  drop  below  the  0.12  threshold  to  cause  an 
increase  in  y.  The  resulting  profile  in  s  is  shown  in  the 
third  panel.  It  is  highly  variable  over  the  profile.  The 
highest  dissipation  is  in  the  lower  altitudes  from  the 
ground  to  twice  the  altitude  of  Mt.  Washington. 

An  expanded  view  of  the  lowest  4km  of  the  atmosphere 
is  shown  in  Figure  6.  There  is  a  region  of  generally 
high  turbulence  in  the  lower  altitudes  with  the  highest 
peak,  actually  off  the  graph  at  1 .25,  there  is  a  minimum 
right  above  1 .4km.  There  is  then  a  region  of  turbulence 
between  1 .6  and  2km  bounded  by  two  peaks,  followed 
by  some  activity,  then  a  drop  to  a  local  minimum  at 
2.5km  followed  by  a  broad  region  of  activity  up  to 
about  3.8km,  with  two  spikes  before  dropping  down 
immediately  above  4km. 

The  next  5km  segment  is  shown  in  Figure  7.  This 
region  is  the  lowest  overall  turbulence  in  the  profile. 

Yet  we  see  several  thin  regions  of  high  dissipation. 
These  are  centered  at:  5.2km,  5.7km,  6.4km,  7.2km, 
7.7km8.75km,  9.3km,  and  10km.  While  many  of  the 
peaks  occur  in  nominally  stable  regions  (Ri  >  0.25),  the 
strong,  narrow  spike  at  7.2km  occurs  in  the  location  of 
minimum  Ri,  and  nearly  zero  potential  temperature 
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gradient  ( --  0).  This  region  at  least  suggests  the 
presence  of  a  nearly  periodic  mountain  wave  structure. 

The  next  5km  segment,  from  10km  to  15km,  is  shown 
in  Figure  8.  This  is  the  tropopause  region,  where  the 
temperature  gradient  is  increasing.  This  causes  higher 
and  generally  higher  values  of  Ri.  In  spite  of  the 
presumed  higher  stability,  the  number  of  regions  of 
high  dissipation  increases,  and  appears  more  random 
than  the  lower  5km  segment. 

The  upper  15km  segment  of  the  sounding  is  shown  in 
Figure  9.  The  turbulence  has  settled  into  a  stratospheric 
pattern.  The  increased  stability,  as  demonstrated  by  the 
higher  and  Ri,  results  in  thinner  layers.  There  are  a 
few  regions  of  sustained  high  turbulence:  the  regions 
between  16  and  17km,  between  23  and  24km,  and 
thinner  regions  centered  on  24.5,  26,  27,  28,  and  29km. 

COMPARISON  TO  LIDAR: 

The  dissipation  rates  estimated  from  the  GroundWinds 
lidar*^  are  compared  to  those  from  thermosonde  launch 
NH_4  in  Figure  10.  The  results  are  generally  in  fair 
agreement,  with  the  notable  exception  of  the  largest 
spike  in  the  thermosonde  estimate  at  1.2km.  Referring 
back  to  the  Figure  6,  this  spike  occurred  at  a  place 
where  the  Brunt-Vaisala  frequency  went  to  zero;  the 


been  in  a  small  region  of  convectively  unstable  air  that 
could  have  been  outside  the  lidar  beam. 

If  we  ignore  that  region  of  disagreement,  the  results  of 
the  two  methods  are  in  fair  agreement.  There  are  two 
minimums  in  the  lidar  estimates  on  either  side  of  3km 
that  are  mirrored  by  the  thermosonde  estimates, 
although  the  lidar  magnitudes  are  lower.  The  regions 
from  about  1.6km  to  2.5,  and  from  3.4  to  4.0  match 
very  well. 

CONCLUSIONS 

Concurrent  velocity  sensing  lidar  and  thermosonde 
measurements  were  made  ihe  evening  of  31  May  2001 
in  the  lee  of  Mt.  Washington,  New  Hampshire.  Data 
from  each  instrument  was  used  to  estimate  the  velocity 
turbulence  dissipation  rate.  The  thermosonde  results 
were  estimated  using  an  extension  of  the  radar  power 
method  described  above.  The  results  showed  very  high 
optical  turbulence  levels  from  ground  level  up  to  4km, 
which  is  twice  the  altitude  of  Mt.  Washington  Summit. 
Above  4km  there  is  a  dramatic  drop  in  optical 
turbulence  followed  by  generally  low  levels  imtil 
activity  increases  near  the  tropopause  and  in  the  higher 
shear  regions  of  the  jet  stream. 

At  publication  date,  the  available  estimates  of 
turbulence  from  the  lidar  are  only  the  first  4km  in 
altitude.  The  comparison  of  the  two  techniques  is  fair, 
with  an  obvious  discrepancy  near  1km,  which  may  be 
due  to  an  isolated  convection  cell,  or  other  horizontally 
inhomogeneous  feature. 
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Figure  10.  Comparison  of  dissipation  rates  as 
estimated  from  GroundWinds  Lidar  data  (dotted) 
and  from  the  Thermosonde  data  (solid). 

Richardson  number  was  near  zero.  Further,  this  spike 
is  in  the  location  of  the  high  Cn  values  associated  with 
capping  inversions  mentioned  earlier.  This  may  have 
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Figure  2,  Profiles  from  thermosonde  flight  Bartlett,  NH  on  31  May  2001.  The  left  panel  shows  air 

temperature  (solid)  and  relative  humidity  (dashed),  the  center  shows  Cn^  along  with  the  smoother  Clearl 
model  profile,  and  the  right  panel  shows  wind  speed  (solid)  and  wind  direction  (dashed). 
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Figure  3.  Profiles  from  the  lowest  5km  of  thermosonde  flight  NH_4.  See  Figure  2  for  variables. 
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Figure  5.  For  thermosonde  launch  NH_4,  profiles  of  the  square  of  the  Brunt  Vaisala  Frequency,  the  gradient 
Richardson  Number,  and  the  turbulent  energy  dissipation. 


Figure  6.  Expended  views  of  Figure  5,  for  the  segment  from  0  to  5km. 


Figure  7.  Segment  from  5km  to  10km. 
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